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Чисельне моделювання процесу кристаллзаци металу 


Рассматривается задача тепловой обработки металла. Строится приближенное решение задачи. Получены 
оценки для температуры и теплового потока. 
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Розглядаеться задача теплово! обробки металу. Будуеться наближений розв’язок задач1. Доведен! ощнки 
на температуру 1 тепловий пот. 
Ключов! слова: процес кристал1заци, математична модель, функця, р1дка фаза, конвекц1я. 


ТЬе ргоШет оЁ фегта! ргосеззте оЁ теа| 1$ сопз14егед. ТБе арргохипайе зоайоп 1$ сопзгасеа. 
Езитаноп$ Юг а 1етрегайиге ап4 еггла] згеатл аге 20. 
Кеу могаб: ргосе5$ оЁ сгуаШтайоп, тафетайса[ то4е|, РапсНоп, Фе Паша рВазе сопуесйоп. 


Распространение тепла в различных средах оказывает большое влияние на ха- 
рактер протекания многих важных для практики процессов. Среди задач, связанных с 
распространением тепла, выделяется класс задач, в которых исследуемое вещество 
переходит из одной фазы в другую с выделением или поглощением тепла. 

Целью данной работы является моделирование процесса кристаллизации 
металла, изучение процесса завершения получения слитка в кристаллизаторе путем 
его вытягивания. 

Рассматривается задача управления информационными процессами при авто- 
матизации технологий тепловой обработки металла, на основе математического мо- 
делирования, анализа статистических данных и теплофизических экспериментальных 
измерений. В качестве источника информации исследуется математическая модель, 
основанная на пространственной задаче Стефана, с учетом конвективного движения 
и примесей в жидкой фазе. 


Постановка задачи 


Пусть ДР = (-1<х<1у<0) — полуполоса, заполненная твердым металлом. 
Обозначим через и(х, у) температуру этого металла. Требуется определить темпера- 
туру и(х, у) по следующим условиям: 


и, чи, аи, =0,(х,у) ЕР, (1) 


УУ 


155М 1561-5359 «Штучний 1нтелект» 2014 №1 153 


Сыпко И.А. 


3С 


и, аи =0,х = 1-9 <у<0, (2) 
и(х,—оо) = 0, (3) 
и, (х,0) =у(х),-1<х 51. (4) 


Здесь @ и @, - постоянные, соответственно, число Пекле и Нуссельта. Решение 
задачи (1) — (4) имеет вид 
Ч со$й„хе“”” 
и(х, у) =, 5 |*(©)соз,645, 
= с0$° Я, °) 
И, (1 т 2 22 ) 


п 


(5) 


2 


@ |@ 
где: = ве + та +1. , п=12,3...А, — положительные корни уравнения А = @, с. 


Отождествим теперь температуру и(х, у) с температурой твердого слитка, на- 
ходящегося в кристаллизаторе при электрошлаковом переплаве. Для вытягивания 
слитка из кристаллизатора поверхность слитка предварительно обогревается тремя 
электронными лучами Й/,Й’, и И’, причем мощность И’, одного из них равномерно 
распределена в центральной зоне ты <х2>ру= 0} ‚ а два других сконцентрированы 
по краям х = +1 [1]. Независимо от того, в каком отношении находится температура 
поверхности слитка с критической температурой Т”, при которой поверхность слитка 
отделяется от стенок кристаллизатора, теплообмен слитка с кристаллизатором осущест- 
вляется по формуле (2). Для получения температуры слитка достаточно положить в 
формуле (5) ъ(х) = (И ,И,,И’,). 

Далее, введем в рассмотрение функционал 


И) = [ма у)-Т) 4 (6) 


Рассматривается задача. Требуется определить поток уъ(х) из допустимого 
множества И , доставляющий наименьшее значение функционалу /(у). Минимизирующая 
последовательность у, строится по формуле у) =у,+=,(у,-у,), параметр =, 
выбирается из условия ши /(у, +5,(у,,-у,)), 0<&=,<1 [2]. В качестве области 
определения функции И берется множество кусочно-постоянных ступенчатых функций: 

УЕ, ,Х, < ХХХ, 1, = СОП5Ё, К = 0},2,...,т. 


При этом формула (5) примет вид: 


а Г(у) = Ци, уу, ,..„Ут 
При численной реализации задачи необходимо учесть ограничение 2500 < (х) < 5000, 
здесь у(х) — мощность потока в единицах МВт/ М. атакже @ = 2,66, @ = 3,05. 
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Способы решения задачи 
Нулевое приближение. Найдем минимум функционала (6), в случае когда 


ио(х, у) = (у), 
с05 ух Ау ель. 


т (ху) =2 те 
Ло Мо (1+ @% 12 °) 


(0 


Минимум функционала (6) находим из условия 


=. = ] (УТУ 0. 


Непосредственные вычисления показывают, что 
2 
с0$° Жо 
№ и 0 (1 + у 2 ) 
Ло 


* 
ый 


у, =4Г - 
эт 2% 


ш24,(-е?®" ; : ш24,@-е^# 
о ЕЕ = о = 
0$ 0$ 
2 (1+ ©, 12 °.) Лоо (1+ 6% 712 °) 
0 0 
Первое приближение. Найдем теперь минимум функционала (6), в случае 


хе’ зтА 

когда и, (х, у) = (№ (ху) + Лоьу)у „где д(ху)=2— 90° ТЯ. 
с0$° д. А, 

и, (1+ @% 12 ) 


1 


Поступая, аналогично тому, как это было сделано в случае нулевого прибли- 


т2, 1-е“) т 24, (1-е) 


жения, получим у =2Т” - - 
о с0$” Жо 2 с0$° Я, 
Лоо (1+ @% 22 ) Ли (+4, 2 
0 1 
где 
р эт? (1-е) эт? 24, 1-е?" +2 812%, 912%, (1-е )(1-е^\) 
_ с05? А. с05? А с05? со5? д 
Ло; (1+ 2 °)? 21 (1+ в, 72 Г)? Лоо + в 12 ^) К, + 2 ) 
1 0 т 


(и 


Далее, имеет место следующая формула: 
2 2иоН 
2% 1- Ио 
Се чи =: 
рр У Зина тт 
1 


0 


т’ 24. 1-е"). 


22 124, 5124 @- е®Н )(1—е®") 
: 1 2 о 
0$ Д 0$“ А, 
м, (1+ @% 2 Г) оо (1+ © 2 °) 


1 0 


| ш 24, (1-е^" ш 2%, (1-е^" . 
2т“, зш 24, (1-е - ‚ 51 „1-е - на”)? 
05 с05 
Др, (1+ @% 12 Г) Лоо 1+ @, 12 °) 
1 0 
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Приближение любого порядка. Аналогичным образом можно исследовать мини- 
мум функционала /(у,) ‚ когда 
п. 0$ 1, Хс0$ АД, (1 е\*" 
е*? +2, >. [ к ( - ) 
сз” Д, А 05° Д, 
Ао (1+ @, 7 °) М, | 1+ 6. р: 


0 А 


у, эт Ду с0$ Их -е^") 


ИкУу 


и„(х, у) =2 


Оценить погрешность предлагаемого метода вычисления минимума функцио- 
нала (6) можно, используя следующее утверждение. 


При достаточно малых значениях @ и при (х, у) е О справедлива оценка 


@,Н оо 
с0$,хс0$ Я, (1-е"*”) емку| < — 


с05’ Я, ен Е 
| 
К 


два 


При доказательстве этого утверждения воспользоваться соотношением: д, =ил+Е,› 
где =, > 0 для всех п[2]. Справедливо также утверждение. 
Пусть выполнены условия ©, > ©2102 ‚ О<@<А, 0<А< =. ‚ О<и, <у(х) <у, 
при хЕ[-1], где у, и у, - некоторые постоянные. Тогда решение краевой задачи 
(1) - (5) удовлетворяет следующим условиям при (х,у) ЕД: 

и, (х, у) < Соехр(иьу) < Сфехр(@у), 
Со ехр(иьу) < и(х, у) < С, ехр(иьу) < С, ехр(@у), 
6+ А(1 + 60/4) ас ЗА 4)1+соз? /А)соз? Ул. 
31+ соз\/ А) 6+ А(1+ соз? /А) 


Для утверждения необходимо сравнить с помощью принципа максимума функции 
и, (х,у) и у(х, у) ‚ где у(х, у) — решение задачи (1) — (5) в предположении, что ъ, (х,0) = у, 


где С, = 


при хе [-1,]. Далее, рассматривается функция /(х, у) =у,(х,у)-и,(х, у), (х, у) Е р 
и доказывается, что /(х, у) >0в Р. Действительно, функция /(х, у) не может при- 


нимать наименьшее отрицательное значение внутри О в силу принципа максимума. 
На вертикальных частях границы х = +1 функция /(х, у) также не может принимать 


отрицательный минимум. В такой точке имели бы Х,(х, у) < 0, между тем Л. (у) = 
=-@/(х,у)>0, х=+1, так как Г(х, у) < 0, по предположению. На бесконечности 
функция Х(х,у) исчезает, т.е. Г(х,-—)=0. На границе у=0, -1<х<1 имеем 
(0) =м, (х,0) -и, (х,0) = у, -у(х) > 0. Следовательно, всюду в р справедливо не- 
равенство и,(х,у)<у,(х,у) при (х,у)е р. Отсюда с помощью интегрирования по 
переменной у следует оценка для функции и(х, у) сверху. Аналогичным образом, можно 
получить оценку на производную и,(х,у) сверху при (х,у)е р. Полученные оценки 
позволяют оценить температуру и(х, у) и тепловой поток внутри области ДР не прибегая 
к решению задачи (1) — (5) [3-6]. 
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Проделанные численные результаты задачи представлены в табл. 1 — 4. 


Таблица 1 — Численные результаты при различных значениях параметров 


3С 


60 № № Т Н о о Г (90) (91) 

1,5 0,9882 | 3,5422 0,9 -3,0 2,783 4,367 1,011 21,997 
1,5 0,9882 | 3,5422 0,9 -4,0 2,783 5,730 5,428 38,576 
1,5 0,9882 | 3,5422 0,9 -6,0 2,783 6,588 -3,038 | 62,490 
1,5 0,9882 | 3,5422 0,95 -3,0 2,937 4,610 1,126 24,509 
1,5 0,9882 | 3,5422 0,95 -4,0 2,937 5,730 -0,559 | 42,982 
1,5 0,9882 | 3,5422 0,95 -6,0 2,937 6,954 -3,385 | 69,627 
1,5 0,9882 | 3,5422 0,9 -3,0 4,5 4,5 16,728 | 23,709 
1,5 0,9882 | 3,5422 0,9 -4,0 4,5 4,5 15,256 | 24,355 
1,5 0,9882 | 3,5422 0,9 -6,0 4,5 4,5 12,397 | 25,275 
1,5 0,9882 | 3,5422 0,95 -3,0 4,5 4,5 15,537 | 23,062 
1,5 0,9882 | 3,5422 0,95 -4,0 4,5 4,5 13,866 | 23,695 
1,5 0,9882 | 3,5422 0,95 -6,0 4,5 4,5 10,718 | 24,696 
1,5 0,9882 | 3,5422 0,9 -3,0 111 11/1 173,683 | 192,101 
1,5 0,9882 | 3,5422 0,9 -4,0 11,1 11,1 175,717 | 198,788 
1,5 0,9882 | 3,5422 0,9 -6,0 11.1 111 173,351 | 202,940 
1,5 0,9882 | 3,5422 0,95 -3,0 11,1 11,1 171,150 | 190,098 
1,5 0,9882 | 3,5422 0,95 -4,0 111 111 172,831 | 196,617 
1,5 0,9882 | 3,5422 0,95 -6,0 11.1 111 170,023 | 200,697 


Таблица 2 — Численные результаты при различных значениях параметров 


60 № № Т Н о о Г (о) (51) 

2,0 1,0769 | 3,6436 0,9 -3,0 3,026 4,949 2,686 26,439 
2,0 1,0769 | 3,6436 0,9 -4,0 3,026 6,212 1,213 47,388 
2,0 1,0769 | 3,6436 0,9 -6,0 3,026 7,642 -1,364 | 78,991 
2,0 1,0769 | 3,6436 0,95 -3,0 3,194 5,224 2,993 29,459 
2,0 1,0769 | 3,6436 0,95 -4,0 3,194 6,557 1,352 52,799 
2,0 1,0769 | 3,6436 0,95 -6,0 3,194 8,066 | -1,520 | 88,011 
2,0 1,0769 | 3,6436 0,9 -3,0 4,5 4,5 16,436 | 20,858 
2,0 1,0769 | 3,6436 0,9 -4,0 4,5 4,5 15,077 | 21,560 
2,0 1,0769 | 3,6436 0,9 -6,0 4,5 4,5 12,270 | 22,543 
2,0 1,0769 | 3,6436 0,95 -3,0 4,5 4,5 15,298 | 20,262 
2,0 1,0769 | 3,6436 0,95 -4,0 4,5 4,5 13,741 | 20,952 
2,0 1,0769 | 3,6436 | 0,95 -6,0 4,5 4,5 10,643 | 22,014 
20 | 10769 | 3,6436 | 0,9 -3,0 11,1 11,1 | 168,455 | 171,436 
2,0 1,0769 | 3,6436 0,9 -4,0 11,1 11,1 171,102 | 178,385 
2,0 1,0769 | 3,6436 0,9 -6,0 11,1 11,1 169,172 | 183,045 
2,0 1,0769 | 3,6436 0,95 -3,0 11,1 11,1 166,052 | 169,559 
2,0 1,0769 | 3,6436 0,95 -4,0 11,1 11,1 168,350 | 176,343 
2,0 1,0769 | 3,6436 | 0,95 -6,0 11,1 11,1 165,972 | 180,926 
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Таблица 3 — Численные результаты при различных значениях параметров 


60 № № Т Н о о Г (00) (91) 
3,0 | 1,1925 | 3,8088 0,9 -3,0 3,352 6,516 7,280 | 45,840 
3,0 | 1,1925 | 3,8088 0,9 -4,0 3,352 8,405 6,004 | 87,001 
3,0 | 1,1925 | 3,8088 0,9 -6,0 3,352 | 10,746 3,378 | 157,183 
3,0 | 1,1925 | 3,8088 0,95 -3,0 3,539 6,878 8,112 | 51,075 
3,0 | 1,1925 | 3,8088 0,95 -4,0 3,539 8,871 6,690 | 96,936 
3,0 | 1,1925 | 3,8088 0,95 -6,0 3,539 | 11,343 3,764 | 175,133 
3,0 | 1,1925 | 3,8088 0,9 -3,0 4,5 4,5 | 20,167 | 18,438 
3,0 | 1,1925 | 3,8088 0,9 -4,0 4,5 4,5 | 19,216 | 19,255 
3,0 | 1,1925 | 3,8088 0,9 -6,0 4,5 4,5 | 16,548 | 20,267 
3,0 | 1,1925 | 3,8088 0,95 -3,0 4,5 4,5 | 19,061 17,873 
3,0 | 1,1925 | 3,8088 0,95 -4,0 4,5 4,5 | 17,904 | 18,668 
3,0 | 1,1925 | 3,8088 0,95 -6,0 4,5 4,5 | 14,927 | 19,741 
3,0 | 1,1925 | 3,8088 0,9 -3,0 11,1 11,1 | 188,971 | 154,728 
3,0 | 1,1925 | 3,8088 0,9 -4,0 111 11,1 | 194,686 | 162,962 
3,0 | 1,1925 | 3,8088 0,9 -6,0 11,1 11,1 | 194,771 | 168,964 
3,0 | 1,1925 | 3,8088 0,95 -3,0 11,1 11,1 | 186,651 | 152,926 
3,0 | 1,1925 | 3,8088 0,95 -4,0 11,1 11,1 | 191,994 | 160,973 
3,0 | 1,1925 | 3,8088 0,95 -6,0 11,1 11,1 | 191,588 | 166,853 


Таблица 4 — Численные результаты при различных значениях параметров 


620 № № Т Н о о Г (5о) (о) 
4,0 | 1,2646 | 3,9352 0,9 -3,0 3,560 8,457 13,134 80,696 
40| 1.2646 | 3.9352 09| -40| 3,560 | 11217 | 12291 | 162,631 
4,0 | 1,2646 | 3,9352 0,9 -6,0 3,560 | 14,954 9,811 | 321,551 
4,0 | 1,2646 | 3,9352 0,95 -3,0 3,758 8,927 14,634 89,911 
4,0 | 1,2646 | 3,9352 0,95 -4,0 3,758 | 11,840 13,694 | 181,203 
4,0 | 1,2646 | 3,9352 0,95 -6,0 3,758 | 15,785 10,932 | 358,271 
40| 1.2646 | 3,9352 09| 30 4,5 45| 26374 | 17.467 
40| 1.2646 | 3,9352 09| 40 4,5 45| 26115 | 18,407 
4,0 | 1,2646 | 3,9352 0,9 -6,0 4,5 4,5 23,830 19,422 
4,0 | 1,2646 | 3,9352 0,95 -3,0 4,5 4,5 25,268 16,899 
4,0 | 1,2646 | 3,9352 0,95 -4,0 4,5 4,5 24,790 17,804 
4,0 | 1,2646 | 3,9352 0,95 -6,0 4,5 4,5 22,171 18,855 
40| 1.2646 | 3,9352 09| 30| 11 11 226/770 | 148,999 
4,0 | 1,2646 | 3,9352 0,9 -4,0 11,1 11,1 | 237,580 | 158,868 
40| 1.2646 | 3,9352 09| -60| 11 11| 241.625 | 166,516 
40| 1,2646 | 3,9352 | 095] 30| ии. 11| 224449 | 147191 
4,0 | 1,2646 | 3,9352 0,95 -4,0 11,1 11,1 | 234,854 | 156,838 
4,0 | 1,2646 | 3,9352 0,95 -6,0 11,1 11,1 | 238,346 | 164,304 
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КЕ5ОМЕ 
ТА. бурКо 
№ итетса Моаейп® орте СгуяаШганоп Ргосеу5 ор Меш 

ВасКогоип4д. Неа ргорагайоп ш уапои$ теда Баз а отеаё шИчепсе оп фе 
сопгзе оР тапу ипропап{ 10 ргасйсе ргосеззез. Атопе Ше ргоШепа$ аззосла ед уу Ше 
ФзтБаноп о Пеаф, Шеге 15 а с1а5$ оРргоетз ш \умеВ Ше шуезижеа зиб$апсе раззез 
Кот опе рВазе 1ю апофег, \иЬ Фе гёвазе ог абзогрНоп оЁВеай. 

Маема 5 ап@ шейод$. Тве ригрозе оЁ 115$ \мотК 15 №ю забуапиме Фе 
тафетайса! то4е| оЁ 1е ргосез$ ог сгубаШтайоп ог те{а|. Тре ргоет оЁ сопо] оЁ 
{есВпо]о21са1 ргосезз оР пеёа[5 {Вегта| ргосеззтс 15 сопз14еге4. Аз шЮгтайоп гезоигсе 
Фе @гее Аилепз1оп сопуесйоп ЗеЁп ргоет ш Па рБазе 15 шуезИеафе4. 

ВезшЕ. ТЬ1$ рарег ех{еп4 10 Ите-4ереп4епе сазе зоте гезиЁ оашеа Бу Ше 
аиог Юг $еа4ду-зае ‘еп рго ет \/ сопуесНоп. То деегтте Фе орйта!| Вега] 
сопай1оп$ оЁ Югтайоп оЁ шеой са]сиЙайоп$ \уеге саге ош ш Фе Нате\мотК оЁ Фе 
тафетайса! то4е| оЁ Фегтла| ргосеззез ш а суппанса| 120%, адаре Юг Фе сазе оЁ 
ВоПо\у 1501. ш Фе то@е| оЁ Паш те{а1 15 роиге4 шю Ше тои! рогйоп$, ап4 а \е4эе 
оЁИ 15 риЙеа о реподсаПу. 

Сопсш$10п. А тафетайса| ое] ус аПо\у$ ю д@егтте 1е егта| то4е ш 
УсВ | 15 дейпед. Зипож{еЯ ргосез$ оЁ сгубаШтаноп ог теёа|, раззше ш зресла| 
теаШогоу, пате[у зе Фе ргосез$ оР сотр!еноп оЁ Ше тесетрЕ оЁ ап шоо{ ш фе 
тош А Бу зтесН1шо. 

Тре ргоет оЁ Шегта| ргосеззте оР шеа| 1$ сопз1Аегед. Тве арргохитае 
зоаНоп 1$ сопзнгисе4. Езитанопз Юг а {етрегашге ап4 Фегтла1 згеат аге 201. 

ТБе г@еуапсе оЁ Пе ргеземе4 \отК 15 4ае Бо 1ю Ше ргасИса1 детап4$ о{ #лтту 
сопего| оЁ Те ргосез$ оЁ сгузбаШтайоп Юг Фе оБес{ ми сотр]ех сеотету. 
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